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Abstract 

We derive a hierarchy of kinetic and macroscopic models from a noisy variant 
of the heuristic behavioral Individual-Based Model of |53j where the pedestrians 
are supposed to have constant speeds. This IBM supposes that the pedestrians 
seek the best compromise between navigation towards their target and collisions 
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avoidance. We first propose a kinetic model for the probability distribution function 
of the pedestrians. Then, we derive fluid models and propose three different closure 
relations. The first two closures assume that the velocity distribution functions 
are either a Dirac delta or a von Mises-Fisher distribution respectively. The third 
closure results from a hydrodynamic limit associated to a Local Thermodynamical 
Equilibrium. We develop an analogy between this equilibrium and Nash equilibia 
in a game theoretic framework. In each case, we discuss the features of the models 
and their suitability for practical use. 
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1 Introduction 

Understanding and predicting crowd behavior is an extremely important issue in our so- 
cieties. Public safety concerns have raised considerably after recent major crowd disasters 
[SI]- Public authorities are challenged by increasingly large crowds attending mass events 
such as sports gatherings. Economic stakes related to crowd management are equally 
high, as increasing the efficiency of pedestrian infrastructures have important returns in 
terms of business. 

To achieve a better comprehension of crowd behavior and increase the reliability of 
predictions, numerical modeling and simulation is playing an ever-growing role. 

A recent review on crowd modeling can be found in [9]. The most widely used crowd 
simulation models are Individual-Based Models (or IBM), such as Rule-Based models 
|61j . mechanical models [351 E3 EH], traffic following models [H], optimal control theory 
models, 00], Cellular- Automata [SS] and Vision-Based models [31 Ell ESI Ell E3 EI] - The 
present paper relies on [53] which is detailed below. Continuum models (CM) are based 
on a fluid dynamics approach [3T)1 13"9~| . Other approaches through optimal control theory 
021 ESI El ES] or exploiting the analogy with car traffic [31 El QUI EH1 EH EH ED] have also 
been developed. For dense crowds, the handling of the volume exclusion constraint has 
led to several specific works [261 123 ED] • Existence theory for some CM of crowds can be 
found in [30] . 

Kinetic Models (KM) are intermediate between IBM and CM. As pointed out in the 
review [9], there are quite few KM of crowds (see an example in [7]). IBM, KM and CM 
constitute a hierarchy of models in the sense that each level can be deduced from the 
previous one by a model reduction methodology. Indeed, KM deal with the one-particle 
probability distributions of IBM. Such a description ignores correlations between the 
particles (which are described by joint probability distributions of k particles for k > 2) 
and is therefore a reduced description of the IBM. CM are deduced from KM by taking 



averages over the velocity variable (see section |4.2. ip . Therefore, CM involve a reduced 
description of the velocity statistics of the KM to a small number of its moments. 

In general, CM or KM are more efficient than IBM for large crowds because their 
computational time does not increase with the number of agents. However, they suffer 
from different drawbacks, such as a reduced validity range due to the necessary recourse to 
closure relations, as detailed below. Nonetheless, CM are invaluable tools for large-scale 
analysis and prediction of crowd behavior. Therefore, it is important to firmly base the 
derivation of CM on their small-scale IBM counterpart. The literature on the derivation 
of CM from microscopic models (IBM or CA) is scarce (see e.g. [U [T51 d3 EI])- The 
present paper addresses this question and intends to propose a hierarchy of KM and CM 
based on the IBM developed in [53] . 

The psychological literature shows that pedestrians can estimate the positions and 
velocities of moving obstacles such as other pedestrians with fairly good accuracy [212] . 
Therefore, the subjects are able to process this information in order to determine the 
dangerousness level of an encounter [66]. Taking these considerations into account, the 
heuristic-based model of [53] proposes that pedestrians follow a heuristic rule composed of 
two phases: a perception phase and a decision-making phase. In the perception phase, the 
subjects make an assessment of the dangerousness of the possible encounters in all the pos- 
sible directions of motion. In the decision-making phase, they turn towards the direction 
which maximizes the distance walked towards their target while avoiding encounters with 
other pedestrians. In this sense, in the game theoretical sense, the pedestrians choose the 
Nash equilibrium as the new direction of motion. Game theoretical approaches of traffic 
have already been considered (see e.g. [62]). 

The goal of this paper is to derive a CM from this heuristic-based IBM. With this 
aim, the time-discrete IBM of [53] is first replaced by a time-continuous IBM and noise 
is added to account for some uncertainty in the pedestrian velocity. From this time- 
continuous IBM, a KM is introduced. The KM describes the evolution of the probability 
distribution function of pedestrians in a phase space composed of position, velocity and 
target direction. For the sake of simplicity, we assume that the pedestrian speed remains 
constant and we discard any slowing down induced by close encounters. We do not develop 
any rigorous theory of the passage from the IBM to the KM [T2l[T6] . 

The passage from the KM to the CM is realized by taking the velocity moments of 
the distribution function. In doing so, some closure relations are needed otherwise the 
hierarchy of moments is infinite. We propose three distinct closure relations. The first one 
assumes a monokinetic distribution function. In other words, the velocity distribution is 
assumed to be a Dirac delta at the mean velocity. Such a monokinetic Ansatz can only be 
valid in the noiseless case but provides an exact solution of the KM. In the second closure 
relation, the velocity distribution function is supposed to be a von Mises-Fisher (VMF) 
distribution [67]. It is adapted to situations where the noise is non-zero. In these first two 
closures, the resulting macroscopic model is a system consisting of a mass conservation 
equation and an evolution equation for the mean velocity of each ensemble of pedestrians 
sharing the same target direction. 

The last closure is based on a formal hydrodynamic limit. It can be performed in 
the restrictive situation where (i) the pedestrian interactions can be approximated by 
spatially local ones and (ii) the interaction region of the subjects is isotropic (i.e. there 



is no blind zone behind the subjects). The closure relies on a Local Thermodynamical 
Equilibrium (LTE) obtained through the solution of a fixed point equation. It expresses 
that each ensemble of pedestrians has found its optimal mean velocity in the midst of the 
other ones, i.e. is a Nash equilibrium in the game-theoretical sense. This example fits in 
the general framework relating kinetic theory and game theory which can be found in [28] 
and which bears analogies with the theory of Mean-Field Games [17]. In a companion 
paper [21], the same methodologies are applied to the model of [56] based on a mechanical 
view of pedestrian encounters. 

The outline of the paper is as follows. In section [2] we review the IBM of [S3]. In 
sections [3] and HJ we successively derive the KM and the CM with the three possible 
closure relations. The resulting models are discussed in |5j Finally, a conclusion is drawn 
in section |6j 

2 The Heuristic-Based model of pedestrian dynamics 

2.1 Principles 

The heuristic-based model of [53] proposes that pedestrians follow a rule composed of two 
phases: a perception phase and a decision-making phase. 

In the perception phase, the key observables are the distance-to-interaction (DTI), the 
time-to-interaction (TTI) and the minimal distance (MD). Let us first examine a binary 
encounter with another pedestrian (see Fig. [1]). 

Definition 2.1 Consider a pedestrian (the subject) at a given location and time moving 
with a given velocity. Suppose that this pedestrian interacts with a single other pedestrian 
(the collision partner) who possibly has a different velocity but whose location at the same 
time is close. In this encounter, we define the following quantities: 

(i) The interaction point is the point on the subject's future or past trajectory where the 
distance to the collision partner is minimal, assuming that both agents move in straight 
line with a constant speed. 

(ii) The Minimal Distance (MD) is this minimal distance between the subject and his 
collision partner. 

(Hi) The Distance To Interaction (DTI) is the distance which separates the subject's cur- 
rent position to the interaction point. The DTI is counted positive if the interaction point 
will be reached in the future of the subject and negative if the interaction point has been 
crossed in the past. 

(iv) The Time-To-Interaction (TTI) is the time needed by the subject to reach this inter- 
action point from his current position (counted positive if this time belongs to the future 
of the subject and negative if it belongs to the past). 

Remark 2.1 (i) If the TTI is negative (i.e. the interaction point has been reached in the 
past and the subject and his collision partner are now moving away from each other), or 
if the MD is above a certain threshold (equal to the subjects ' diameter, possibly augmented 
by some safe-keeping distance), then, no interaction occurs and the DTI is set to infinity. 



(ii) Because the subjects have supposedly perfect knowledge of their own and partner's 
positions and velocities, we assume that they are able to estimate the DTI, TTI and MD 
with perfect accuracy. 

We now define the DTI and the TTI when several collision partners are present (see 
Fig. |2J). We have: 

Definition 2.2 When the subject is interacting with several collision partners at the same 
time, the subject's global DTI is the minimum of the DTI of all binary encounters. We 
denote it by D(w) if w is the velocity of the subject. 

The decision-making phase consists in changing the current cruising direction u to 
a new cruising direction u' . It is the outcome of an optimization process. From the 
knowledge of the DTI in each cruising direction, the subject chooses the direction which 
maximizes the DTI, while keeping his direction of motion close to his target direction. In 
|53j . the decision making phase is performed at discrete times separated by equal time 
intervals At. In this phase, the cruising direction is updated through the following local 
optimization procedure. Without any obstacle, the subject would choose a target direction 
a (a is a unit vector of IR 2 ) and cruise with speed c. Therefore, after a time interval At, in 
the absence of obstacles, he would find himself at position X? = c At a (assuming that the 
origin of the coordinate system is placed at his current position). The point X? is called 
the target point. Now, in the presence of obstacles, the subject cruising in the direction 
w (w being a unit vector in M 2 ) will estimate impossible to move a distance larger that 
the DTI. Therefore, choosing the cruising direction w will place the subject at a position 
Xe(w) = D{w)w. The point X E {w) is the estimated point reached in the direction w. 
The decision making consists in choosing for new cruising direction u' the direction such 
that the estimated point Xe{u') is the closest to the target point Xt, among a set of 
test directions w belonging to the vision cone C u about the subject's current direction of 
motion u. Therefore, u' is determined by 

u' = arg min \X T — X E (w)\, (2.1) 

wec u 

where arg min denotes the point that realizes the minimum. Such a realization of the 
minimum may not be unique, but we will discard this possibility as non-generic. This 
decision-making phase is illustrated on Fig. [3J 

This decision-making rule implicitly states how binary interactions are combined. This 
combination is not a mere superposition, as in the classical social force model [35l 13?] 138] . 
but a highly nonlinear operation involving the solution of an optimization problem. In 
[52] , it is shown that this model provides a better account of some of the most striking 
emergence phenomena in crowds, such as spontaneous lane formation in bidirectional 
motion. 

In the next two sections, we make all these considerations mathematically explicit. 

2.2 Perception phase 

We consider a pedestrian i located at a position £,(£), with a velocity V{. He interacts 
with a collision partner j located at a position Xj(t) which has a velocity Vj. Figured] 
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Figure 1: The Minimal Distance MD and the Distance- To- Interaction (DTI). The MD 
is the smallest distance which separates the two pedestrians i and j supposing that they 
cruise on a straight line at constant velocities t>j and Vj. The point on pedestrian z's tra- 
jectory where the minimal distance is attained is the interaction point Xjint of pedestrian % 
in his interaction with pedestrian j. The MD is the distance between xa Qt and Xj int . The 
DTI is the distance which separates the current pedestrian position x, to the interaction 
point Xjint- The Time-To-Interaction (TTI) is the time needed by pedestrian i to reach 
the interaction point from his current position. Clearly, TTI = DTI/|t>j|. 
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Figure 2: Perception phase: the Distance- To-Interaction (DTI) of a given pedestrian in 
the case of several simultaneous encounters is the minimum of the DTI of the individual 
encounters. In this figure, the DTI of pedestrian i with pedestrian j (denoted by DTI(j)) 
is smaller than that with pedestrian k (denoted by DTI(fc)). Therefore, the DTI of 
pedestrian i is DTI(j). 
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Figure 3: Decision-making phase. The new cruising direction u' is chosen such that the 
estimated point Xe{v!) in this direction is the closest to the target point Xt- The tested 
cruising directions w (with associated points Xe(w)) are restricted to the vision cone C u 
of the pedestrian (where u is the current cruising direction). 



gives a schematic picture of the geometry of the collision. The goal of this section is to 
compute T int , -Di nt , -D min , respectively the TTI, DTI and MD of walker i in his interaction 
with pedestrian j (see Definition 12.11) . 

Lemma 2.3 We have: 



Tint 
Dint 

J-/ rn.ir). 



[Xj - Xj) ■ [Vj - Vj) 

\Vj -Vi\ 2 

(Xj - Xj) ■ (Vj - llj) 

\Vj -Vi\ 2 



V; 



Vj-Vi ,2\ 1/2 



Xj X{\ \.\*^i *^i) 



\ Vj - Vi 



(2.2) 
(2.3) 
(2.4) 



Proof: The distance D(t) between the two particles at time t is given by 
D 2 (t) = \xj + Vjt - (xi + Vit)\ 2 



\v j -v i \ 2 (t + 



(Xj -Xi)- (Vj - Vi) 



~\~ \X j X{ 



((Xj-Xi) '(Vj-Vi))' 



\Vj-Vi\ l J \Vj-Vil 2 



(2.5) 



denoting by Xi and Xj the positions of the two particles at time 0. This quadratic function 
of time is minimal at the time t = r- mt given by (J2.2J) . which gives the value of the TTI. 



MD =\. 

this distance 



TTI = this distance 
divided by \vj — u,-| 



Figure 4: Geometry of a collision: The TTI is the projected distance of the two pedestrians 
on the direction of the relative velocity Vj—Vi, divided by the norm of this relative velocity 
\vj —Vi | . The DTI is the TTI times the velocity of the pedestrian. The MD is the projected 
distance of the two pedestrians on the normal direction to the relative velocity. 



Then, the DTI D int of particle i is obviously given by the distance traveled by this particle 
during the TTI, i.e. D int = r int \vi\. This leads to (j2.3p . Finally, the MD D min is given by 
the minimal value of ( 12.5(1 . i.e. D min = D(r int ), which leads to ( 12.41) . ■ 

If ( 12. 2 p and (12. 3p give negative values for the TTI and DTI, it means that there is no 
threat of collision in the future times, as pedestrians are walking away from each other. 
Therefore, some interaction occurs in the future if and only if (xj — Xj) • (vj — V{) < 0. 
Furthermore, if the MD is larger than a certain threshold R identified as the diameter 
of the individuals, plus a certain safe-keeping distance, the interaction will no longer be 
perceived as a collision threat. In both cases, the DTI and TTI are set to infinity. With 
these additional features, we now define the DTI and TTI as 

Definition 2.4 We define: D int and r^t as: 

T int- \vj-vtf ' U int- \vj-vtf l V *l> I ^6) 

if (\xj - Xi\ 2 - ((xj - Xi) ■ |^J|) 2 J < R, and (xj - x { ) ■ (vj - v$ < 0, 
Ti nt = +00, D int = +00, otherwise . (2.7) 



2.3 Decision-making phase 

The collision avoidance model of [S3] uses the elements of collision perception reviewed in 
section [2721 In this model, the decision of a new cruising direction taken by the pedestrian 
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reflects the balance between two antagonist goals: collision avoidance on the one hand 
and maintenance of the target direction on the other hand. The goal of the present paper 
is to investigate the role of the cruising direction. Consequently, we discard the variations 
of the cruising speed. We assume that all pedestrians move with constant speed equal to 
c. Therefore, \vi\ = \vt\ = c and we let 

Vi = CUi, Vj = CUj, \Ui\ = \Uj\ = 1. 

This assumption prevents us to take into account one of the features of the model of 
[5*3] . namely that pedestrians slow down or stop in case of very close encounters. This 
restricts the validity of the present paper to low densities. In the present paper, we will 
also discard fixed obstacles. 

The model follows the lines of [53] , with some simplifications of the expressions of the 
collision avoidance response. We assume a time discrete model with time steps At. During 
this time step, the pedestrian moves according to the vector ViAt. Then, he updates his 
velocity and adopts a new velocity. For this purpose, he explores all possible velocity 
directions w and computes the minimum of the DTI with the other pedestrians in the 
direction w. Let us denote by Dij{w) the DTI with pedestrian j assuming that i moves 
in the direction w. If there are no close encounters, we let this quantity be equal to the 
distance traveled by the pedestrian during At, i.e. cAt. In any case, we limit D^w) by 
this quantity. Then, according to (j2.6p . (12 .7p . we have: 



/ 



Dij(w) 



(Xj — Xi) ■ (Uj — w)\ \ 

min ( J j — — L , cAt j if (Xj — Xj) • [Uj — w) < 0, 



\Uj — w\ 



and ^-^-((z,-*,)-^) 2 ^ 2 , ( 2 - 8 ) 
c At otherwise. 



For physical consistency, we should have R < cAt, as the typical diameter of a pedestrian 
should be much less than the distance traveled by a subject between two velocity updates. 
Then, we define the minimum Di(w) of all the DTI by taking the minimum of Dij(w) 
over all partner pedestrians. The anisotropy of human vision is taken into account by 
restricting the set of partner pedestrians to those belonging to the vision cone of pedestrian 
i. Introducing a threshold number k G [0, 1], this vision cone is centered at Xi and has half 
angle cos -1 k about the direction itj. The minimal DTI of the ith-pedestrian is therefore 
defined by : 



*>-/ -1 fJj i 



AH=min{AiH, j = l,...,N, j^i, -^ l -7-Ui>K}, (2.9) 



*Aj q •-*-"). 



where N is the total number of pedestrians. Finally, the new direction of motion u' { of the 
z-th pedestrian is the direction w that minimizes the distance between the point reached 
after traveling a distance Di(w) in this direction and the point reached after traveling a 
distance cAt in the target direction aj. Therefore, the new direction of motion u[ is found 
by solving the following minimization problem: 

u[ = arg min \Di{w)w — c Ata,j| 2 , (2-10) 

MISS 1 | W-Ui>K 



where again, the test directions w are restricted to the vision cone of pedestrian i. We 
denote by S 1 the set of vectors of M 2 of unit norm. 

Remark 2.2 We note that the minimization problem $2.10\) is not convex and may have 
multiple solutions. We will discard this occurrence as non-generic. 

2.4 Summary of the Heuristic-Based IBM model 

We now consider a collection of N point particles with positions x™ G IR 2 , velocity direc- 
tions u™ G S 1 at time t n = nAt and target direction a, G S 1 , and i G {1, . . . , N}. The 
dynamics is as follows: 

*r +1 =*? + cAt<, (2.11) 

< +1 = arg min \D^ +1 {w)w - c Ata^ 2 , (2.12) 



with 



x n+1 — x n+1 



Dr L M = min{D^(w), j = l,...,JV, j? i, , J n+1 _ * n+ll •<>«}, (2.13) 

Xj x i 



and 



< ■ \x n , +1 - x^ 1 ) ■ {u™ - w)\ 



D?+ l (w) 



mm 



|2 



cAt). 



\u-—w\* J 

if (x] +1 - x™ +1 ) ■ (u] -w)<0, 



and |^ +1 - x7 +1 \ 2 ~ ((x] +l - xr 1 ) ■ $%y < R\ 
c At otherwise. 



(2.14) 



We now make some comments on the position update rule (12.111) . Since the pedes- 
trian can only walk a distance D{u™) in the direction u" before colliding with another 
pedestrian, it would appear more sensible to use the formula a;" + D(u™) u™. However, in 
the present model, the pedestrian speed is supposed equal to one. Therefore, this update 
can only provide the position at an intermediate time t n + D(uf)/c < t n+1 . This leads to 
position updates at different times for the different pedestrians. In order to derive a time 
continuous model, it is more convenient to keep position updates at constant time-steps, 
which justifies the choice made in (12. lip . In the limit At — > in (12. lip (but keeping At 
finite in (I2.12p . (I2.14p ). the quantity cAt tends to zero and eventually becomes smaller 
than D{u™). Then, the objection formulated at the beginning of this paragraph disap- 
pears. In the next section, we first propose a time-continuous model based on this limit, 
and then deduce a mean-field kinetic model from it. 

3 Mean-field kinetic model 

3.1 Methodology 

The goal of this section is to propose a time and space continuous kinetic model (KM). 
With this aim, we first convert the previous time-discrete IBM into a time-continuous 
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one. This conversion consists in replacing the sudden change of the velocity every At 
time intervals, by a continuous one. 

The difficulty with writing such a time- continuum model comes from the 'roughness' 
of the rules of the time-discrete IBM. For this reason, we regularize the time-discrete 
dynamics in two ways. First, in the perception phase, we replace the distance to the closest 
encounter by an average distance to the possible encounters in some interaction region. 
We propose the use of an harmonic average which closely approximates the minimum 
used in the original model. The use of averages over certain interaction regions is found 
in many classical swarm models, such as [21 [2U1 123 EH E31 ES], but the introduction 
of harmonic averages is new. Second, in the decision-making phase, the jump to the 
direction of motion which maximizes the distance walked towards the target is replaced 
by a continuous directional change determined by a velocity potential. This supposes 
that the subjects choose their new cruising direction close to the previous one, which 
looks realistic. 

A final modification of the IBM is to add some noise in the pedestrian velocity updates. 
This noise accounts for various effects such as the uncertainties in the estimations of the 
interaction partner velocities, the variability of the subjects' responses to interactions, 
the possibility that the subjects react to some unpredicted stimuli, etc. For KM, the 
introduction of noise in the particle velocity update results in diffusion in velocity space 
which produces solutions with smooth velocity profiles. This has important consequences 
for the derivation of CM from KM, as it supports the use of smooth macroscopic closures. 
Such smooth closures will be at the heart of the VMF closure and of the hydrodynamic 
limit methodologies which will be described in Sections 14.31 and 14.41 respectively. 

3.2 Modified time-continuous IBM 

We consider the following time-continuous stochastic model for the pedestrian positions 
Xi(t) and velocity directions Ui(t): 

^ = cn,(t), (3.1) 

dui = Fi(t)dt + ((V2dodBi(t)) ■ uf)uf, (3.2) 

where F{(t) is a force term and d is the noise intensity (supposed uniform among pedestri- 
ans). The term dBi(t) stands for the standard white noise and the symbol V means that 
the stochastic differential equation must be understood in the Stratonovich sense. The 
force term Fi(t) is constructed below in such a way that it remains orthogonal to u»(t), i.e. 
Fi(t) ■ Ui(t) = and the noise term V2ddBi(t) is projected onto the orthogonal vector uf- 
to Uj. These facts, together with the use of the Stratonovich definition of the Stochastic 
Differential Equation, maintain u, on the one-dimensional unit sphere i.e. |Wi(t)| = 1, 
provided that |w.;(0)| = 1 initially [35]. 

The force term is defined as follows. First, we replace the 'min' in (I2.13P by an average 
over neighboring particles located in the z'-th pedestrian interaction region. We choose an 
harmonic average, which has the property to give large weights to the small values of the 
quantity to be averaged. In this way, the harmonic average mimics closely the outcome 
of the 'min' operation in ( I2.13p . The z-th pedestrian interaction region is defined as the 
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angular sector centered at Xi(t), with axis v,i(t), semi-angle cos l k and radius Si(t). The 
set Si(t) of subjects belonging to the z'-th pedestrian interaction region is: 



S i (t) = {jE{l,...,N} | \x 3 {t)- Xl {t)\<m and ff% ^ • Ui (t) > k}. (3.3) 



Xj(t) -Xj(t) 
\Xj{t) - Xi{t)\ 



The angle cos -1 k, is the semi-angle of the human vision cone (say typically it/ 2, i.e. k — 0). 
The value of <5j(£) is linked to the local inter-particle distance and will be estimated later 
on. The number of elements of Si(t) is denoted by #Si(t). 

We then consider the harmonic average of the elementary DTI with all collision part- 
ners in the interaction region: 

Dr 1 {w,t)=max{—^ m £ D^(w,t),±}, (3.4) 

where the 'max' has been introduced to bound the average for reasons that will be ex- 
plained below and where the DTI is defined like in ( 12.141) : 



/ 



D^(w,t) 



Uj(t) — w 
mm ' 



2 



\(xj(t) -Xi(t)) ■ (uj(t) -w)\ 
if (xj(t) - Xi(t)) ■ (uj(t) - w) < 0, ( 3 5 ) 

and \ Xj (t) - x t (t)\ 2 ~ {(xj(t) - Xi(t)) ■ ^g^) 2 < R\ 



otherwise. 



The quantity £ is a lower cut-off for Dij because the elementary DTI can be arbitrarily 
small. In reality, if the DTI is too small, the pedestrian lowers his velocity or even 
stops. This feature is not taken into account in the present model, where we only allow 
directional changes. Therefore, to cope with this situation, our model pedestrian would 
have to develop very large angular accelerations, which is unrealistic. The parameter £ is 
introduced to bound the forces and thus prevent the dynamics to become too singular in 
this situation. In the situations where the elementary DTI are large (which corresponds 
to the second alternative of (13. 5p ). we just set them equal to oo, so that they are not 
taken into account in the average ( 13. 4p . which computes the global DTI. The bound of 
the global DTI by cAt is realized by the parameter L as described now. 

Indeed, the quantity L stands for the distance walked by the pedestrian between two 
velocity updates (i.e. L = cAt in the discrete model). Of course, if there are no collision 
partners (i.e. Si(t) = is the vacuum set), or if the elementary DTI with the available 
collision partners are large, the pedestrian will be able to walk this distance L without 
performing a velocity update. Therefore, we bound Di by L thanks to the 'max' in (I3.4J1 . 
In practice, we have £ ~ R ^ L. Indeed, the lower cut-off for the elementary DTI is of 
the order of the size of a subject, while the free walking distance between two velocity 
updates is much larger. 

We now define the i-th pedestrian potential function $*(w) for unit vectors w G S 1 
by: 

$ t i (w) = -\D l (w,t)w-La i \ 2 . (3.6) 
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The coefficient (kL 2 )~ l gives the order of magnitude of the potential and of the force. By 
(13. 2 p and the fact that the velocity u is dimensionless, the force and consequently (kL 2 )^ 1 
have the physical dimension of a reaction rate. Therefore, we can view the quantity 
(fcL 2 ) -1 as providing the typical magnitude of the pedestrian reaction. The force Fi(t) is 
defined by 

Fi(t) = -V^Kmit)) (3.7) 

= -k ( (A(ui(*),t) -Loi(t)-Ui(t)) V w Di(ui(t),t) 

-Di(ui(t),t) L (oi(t) ■ w i (t) ± )w i (t) ± ) . (3.8) 



We note that gradients of functions defined on S 1 are tangent fields to S 1 . Therefore, 
by formula (13. 7J , Fi(t) is orthogonal to Ui(t) as it should. This is reflected in (I3.8p . The 
first term is proportional to V w Di(ui(t),t), while the second one is proportional to uj-, 
and both are orthogonal to Ui(t). 

Definition (13. 7p reflects the fact that, under the force F i: the pedestrian decreases his 
potential $*. Therefore, the pedestrian turns towards the direction of the local minimum 
of the attraction basin of $* to which he belongs at time t. This rule can be viewed 
as a local version of the global minimum rule set up by (I2.12p . Of course, this local 
minimum may not be the global one expressed by (12.121) . However, whether in actual 
life, a pedestrian spontaneously chooses the global minimum or a local one close to his 
current direction of motion is not clear. Therefore, to our opinion, this local rule is as 
legitimate as the global one, until experiments can clarify this point. Of course, the two 
rules coincide if the local minimum is equal to the global one. So the discrepancy between 
them may be quite small in practice. 

3.3 Mean-field kinetic model 

We now introduce a statistical description of the system. Instead of using the 'exact' 
positions, velocities and preferred directions of pedestrians, we rather describe the system 
in terms of the probability distribution f(x,u,a,t). Specifically, f(x,u,a,t)dxduda is 
the probability of finding pedestrians in a small physical volume dx about point x, within 
an angular neighborhood du of velocity direction u, and within an angular neighborhood 
da of preferred direction a at time t. We recall that x £ 1R 2 , u, a £ S 1 . If there are 
no interactions between the pedestrians, i.e. if the acceleration term F is due to purely 
external causes, f(x, u, a, t) can be rigorously proved to satisfy the following Fokker-Planck 
equation: 

d t f + cu ■ VJ + V u ■ (Ff) = dA u f. (3.9) 

This equation is a consequence of Ito's formula of stochastic calculus. The left hand- 
side is a transport operator. It expresses the material derivative of / in the phase space 
spanned by (x,u), due to the motion of the particles with velocity cu and acceleration 
F. The right-hand side is a velocity diffusion term which comes from the velocity noise. 
Let 9 be the angle between u and the first coordinate direction. Then, u = (cos 6, sin 9), 
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u = (— sin#, cos#), and each term of (I3.9P is written as follows: 

u-VJ = cos9d x J + sin9d x J, V u ■ (Ff) = d e (F e f), A u f = d 2 e f, 

where the force term F = Fqu 1 - is by definition orthogonal to u because \u\ = 1. 

We note that there is no operator acting on the a-dependence of /. This is because 
we assume that the target direction a is a quantity attached to the agents which is not 
changed with time. This assumption could easily be modified to take into account a 
possible change of the target direction with the motion of the pedestrians. However, even 
with this simplifying hypothesis, the statistics of target directions has a definite influence 
on the dynamics through the interaction force described below. 

Here, the acceleration term F is not due to external forces but to interactions between 
the particles. So the rigorous derivation of (I3.9P is more difficult and is left to future work 
(see e.g. |12j). The acceleration F is coupled to / through continuous equivalents of ( 13 IT]) , 
and is written: 

F(x, u, a, t) = -V w $(a. )U ,o,t)(u), (3.10) 

where $(x,u,a,t)(w) is the potential of a pedestrian at time t located at x with velocity u 
and target velocity a. The potential is a function of the test direction w. It is given by 



k 



®{x,u,a,t)(w) = -\D iX)U ^(w)w-La\ 2 , (3.1T 



in terms of the DTI D^ x , u ,t){w) of pedestrians located at position x at time t with velocity 
u in the test direction w. 

To compute the DTI, we first define the interaction region of such a pedestrian by: 

S(x, u, t) = {y e M 2 | \y - x\ < 5{x, t) , p ~ - -u>k}, (3.12) 

where S(x, t) will be estimated later on. Then, the continuous equivalent of ( 13. 5 j) leads to: 

n _x , , iIyeS(x,u,t)I(v,b)eT^-\y-x,v-w)f(y,v,b,t)dydvdb i| 

F) t „. Aw) = max < —^ — 7 . — — — — - — , — > .3.13 

(X ' M ' 4) 1 f y eS( x ,u,t)f {v , b) eT>f(y,v,b,t)dydvdb 'l/ 1 

We have denoted by T 2 the two-dimensional torus T 2 = S 1 x S 1 . In ( 13. 13ft . the quantity 
D(y — x,v — w) is the elementary DTI of a pedestrian located at position x and velocity 
w in the encounter with a particle located at y and having velocity v. It is given by: 



D 1 (y — x,v — w) 



\v — w\ 2 1 

mm ' 



\(y — x) ■ (v — w)\ I 
if (y - x) ■ (v - w) < 0, (3.14) 

and 1^ - ^i 2 - ((2/ - x) - ^E^,) 2 < ^ 2 , 

otherwise. 



The significance of formulas (I3.13P and (I3.14p and the roles of the parameters £ and L are 
the same as in the time-continuous IBM of section | 
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Collecting A3. 13[) and (I3.14p allows us to compute the potential $( x ,u,a,t)( w ) given by 
(13. lip . After computing the gradient, the force (J3.10P has the expression: 



,*) 



F(x,u,a,t) = -k ((D(a. iUjt )(u) -La- u) V W D^ U - 

-D (XiU:t) (u)L(a-u x )u x y (3.15) 

Now, we can provide an estimate of 8(x,t). As the density increases, the mean inter 
particle distance decreases like N~ l l 2 (x,t) where N(x,t) is the local density: 



N(x,t) — I f(x,u,a,t)duda. (3.16) 

J{u,a)eT 2 

Therefore, the DTI should decrease in the same proportion. One way to achieve this 
scaling is by taking S(x, t) ~ N" 1 ^ 2 (x, t). Indeed, since D(y — x,v — w) is of the order of 
\y — x\ (by (13.141) ). the average DTI D^ u ^{w) is of the same order. And since \y — x\ < 5, 
we obtain the expected scaling of D( XjUi £) (w) like N~ 1 ^ 2 (x, t). In practice, we need to take 

5(x,t)=CN- 1/2 (x,t). (3.17) 

with C sufficiently larger than 1 to ensure that the estimate (13.131) will take into account 
enough pedestrians. 

Finally, the KM consists of the kinetic equation (13.91) . with the acceleration computed 
through (13151) . 

3.4 Mean-field kinetic model: discussion 

The mean-field model expresses how the statistical distribution of the pedestrians in 
position, velocity and target direction evolves with time. This evolution combines a 
transport operator (left-hand side of (13. 9p ) which describes pedestrian motion towards 
their target direction and collision avoidance, and a velocity diffusion operator (right- 
hand side of (13. 9p ). which models velocity uncertainty. The pedestrian speed c is supposed 
constant because the model focuses on directional changes only. Directional changes are 
modeled through a force term F (13.151) . which describes how pedestrians find the best 
compromise between their target and the necessity of avoiding pedestrians passing by. 

The force F is tailored to decrease the potential function ^f x ,u t t)( w )- This potential 
describes how well the target point is approached when the pedestrian (initially located 
at position x, velocity u and target velocity a at time t) moves in direction w (formula 
(13.111) ). For a set of test velocities w, the pedestrian computes his DTI D(y — x,v — w) 
with a pedestrian located at y with velocity v (formula (13.141) ) and averages it over all 
pedestrians located in his vision cone S(x,u,t) (formula (13.121) ). giving rise to D^ x>u ^(w) 
(formula (13.131) ). This averaged DTI provides him with an estimate of the distance he 
can move in the direction w and allows him to compute his potential &(x,u,t)(w). Finally, 
the pedestrian turns to ensure the decay of the potential and to get closer to his goal 
(formula (I3.10p ). The interaction term is spatially non-local, through (I3.13p . In the next 
section, we derive a spatially local approximation of this non-local term. 
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3.5 Mean-field kinetic model with local interaction 

If we observe the system at a large distance, the various length scales involved in the 
interaction terms appear to be small. Therefore, under this assumption, it is legitimate 
to assume that there exists a small dimensionless quantity 77 <§C 1 such that 



5 = 7)8, R = TjR, L = r]L, £ = iji, 



(3.18) 



where all 'hat' quantities are assumed to be 0(1). Simultaneously, we assume that the 
pedestrian reaction rate remains 0(1). We recall that the pedestrian reaction rate is 
measured by the coefficient kL 2 (see discussion after Eq. (13. 6p ). This assumption implies 
that 

rfh = k = 0{l). (3.19) 

We introduce the change of variables y = x + 7/f , with £ G M 2 , in (13.131) and keep only 
the leading order terms in the expansion in powers of rj. In this scaling Eqs. (I3.9P and 
(13.101) are unchanged, except that all unknowns f 71 , F v , <& v now depend on rj. Then, the 
condition y G S(x, u, t) is equivalent to the condition £ G C u K $, where 



C 



U,K,8 



{U 



If I < 5 and 



r «>K}. 



We have 



with 



D (y — x, v 



w\ 



V 



w) 



(3.20) 



mm 



D-\^v-w)={ 



\v — w\ 



1 



£>(v-w)\ £ 
otherwise. 



if f • (v - w) < 0, 

^d |£| 2 " (£ ■ ffif < R 



(3.2i; 



Consequently, 



with 



^~i,t)H 



D (Lt)( w ) 



1 



-D 



(x,u,t) 



»» 



hec Jivtof&D 1 (£ 1 v-w)f(x + rj£,v,b,t)d£dvdb 1 

u,M f T 77 - t ...,,, „ , - } ■ (3.22) 

J fee _-J< vm Y*f\ x + rtiV,b,t)dtdvdb 



Finally, we have 



®(x,u,a,t)(w) 



k, " - ,2 

77|- D (x,n,t)(w)w-La| . 



Since we look for a local approximation, we assume that rj <C 1. The distribution / 
is assumed to evolve only on the large scale. Therefore, in the Taylor expansion of (I3.22p 
with respect to 77, we may keep only the leading order term and neglect the higher order 
ones. As a result of this approximation, f disappears from the arguments of the function 
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/ in both the numerator and denominator. The integration with respect to £ can thus be 
performed beforehand, leading to the quantity A~ l Ju, v — w) defined by 



A >"- ) = ssfcj/^ 1 6 - , «--»)* 



(3.23) 



and Area(C u K §) is the two-dimensional area of C u K $. This leads to the following expression 
of D, dropping all the hats for simplicity: 

Again, 5(x,t) is linked to the total density through ( I3.17p . Graphical representations of 
D~ l and A K) g can be found in Figs. |5] and |6] respectively. They illustrate that the function 
A Kj< 5 only depends on u ■ ,^~^ and \v — w\ (i.e. two real variables) while a general function 
of (u,v — w) depends on a vector of S 1 and a vector of K 2 , i.e. three real variables. This 
is due to the fact that D itself only depends on £ • {v — w) and \v — w\. The function A K) s 
can be numerically computed a priori. 

We note that the expression of A K) ^ simplifies in the special case k = —1. In this case, 
there is no blind zone: all their collision partners in the disk 

B S = {£ e M 2 | |f | < 5}, (3.25) 

are taken into account by the pedestrians. Consequently, the averaging (I3.23P is over the 
whole disk B$ (see Fig. [7]). The dependence of A_ 1)( 5 upon u disappears. The resulting 
function, denoted by As(\v — w\), is given by: 

A ^-^ = ASaW/ Bi 6 " 1 «'"- ro »^ (3 - 26) 

Additionally, it is an elementary matter to remark that when \v — w\ — > 0, we have 
A^ 1 (|t; — w\) ~ |v — iy| In ' v ~ w ' . Another consequence of this simplification is that the 
potential $ does not depend on u. It can be simply written $> x ,a,t(w). This simplification 
will be exploited in the hydrodynamic limit (see section |4~4|) . 

We summarize this section: due to the assumption that / evolves on the large scale 
only and thus can be taken constant in the interaction region of a given pedestrian, the 
interaction force only depends on / at that location. This local approximation scaling 
(I3.18p . ( I3.19P leads to the kinetic model (I3.9P with a local evaluation of the force. The 
force is still computed from the potential (13. lip through (I3.10p . However, the evaluation 
of the DTI is now given by a local velocity average ( 13.24J1 . where the velocity convolution 
kernel A~^(w, v — w) can be analytically computed. In the simpler case where there is no 
blind zone, the kernel reduces to a function of \v — w\ only, and the potential ^ x ^ t {w) 
does not depend on u. 
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\v — w\ 

Figure 5: The function D(!;,v — w). For a fixed value of v — w (in red), the parallel and 



transverse components of £ are respectively denoted by £y = £ ■ 



\v— w\ 



and £± = £ — £|| ■ The 



domain of definition of £ — > D(£,v — w) (domain where it is finite) is characterized by 
£|| < and |£j_| < R (see (j3.2ip ) and is the shaded blue area on the figure. The function 

£ — > D(£,v — w) is constant along all vertical segments and has value , J 1 , except in 
the dark blue area where it is constant equal to £. The transition happens along the line 
£|| = — £\v — w\ (the vertical blue line on the figure). 




£u 



Figure 6: The function A Kj g(u, v — w) is obtained by averaging the function £ — > £) -1 (£, v — 
w) (see Fig. [5]) on the cone C Ut s (represented by the green shaded area). The blue-green 
shaded area is the intersection of the cone C UiKt s and the support of £ — > l) -1 (£, v — w). 
It highlights the fact that the function A Kt § depends on u ■ , V Z W \ an d \ v ~ w \- 



The disk B* 




Figure 7: Case k = — 1 (i.e. there is no blind zone: all collision partners in the disk B$ 
of radius 5 are taken into account). The function A$(\v — w\) is obtained by averaging 
the function £ — > Z) _1 (£,t> — w) (see Fig. [5]) on the disk B$ (represented by the green 
shaded area). The blue-green shaded area is the intersection of B$ and the support of 
£ —7- £) _1 (£, v — w). The function As depends only on \v — w\. 



4 Macroscopic model 

4.1 Necessity of a closure Ansatz 

Macroscopic models are obtained by taking averages of functions of the particle velocity 
u over the distribution function f(x,u,a,t). The resulting macroscopic quantities are 
e.g. the density p(x, a, t) or the mean velocity U(x, a, t) of pedestrians at position x with 
target direction a at time t: 

p(x,a,t) = / f(x,u,a,t)du, (4.1) 



1 f 

U(x,a,t) = — r / f(x,u,a,i)udu. (4.2) 

p{x,a,t) J u& i 

It is necessary to keep the dependence of the macroscopic quantities over the target 
direction a. Indeed, in general, the statistics of the target directions is not known or may 
change from situation to situation. In situations where the statistics of target directions 
is known and does not change with time, it is also possible to introduce the total density 
N(x,t) (already met at (13.161) ) and the average velocity V(x, t) of the pedestrians at 
position x and time t, irrespective of their target direction. The latter is defined by: 

V(x,t) = / f(x,u,a,t)ududa. 

N[x,t) J MeJ 2 

In the present work, we will only consider CM which retain the statistics of target direc- 
tions. 
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To pass from the KM (I3.9P to a CM for the quantities p and U, the most direct 
method is the moment method. It consists in integrating the kinetic equation (13. 9p with 
respect to the particle velocity u, after pre-multiplication by polynomial functions of u. 
Unfortunately, in general, this method does not yield a closed model for p and U because 
higher order moments (i.e. integrals of higher order polynomials of u) may be involved 
in the resulting system. These higher moments need to be expressed in terms of p and 
U through a suitable closure relation. Closure relations are usually provided through an 
Ansatz which expresses / itself as a function of p and U. The validity of this Ansatz is 
subject to caution. When dissipative phenomena are present, such as in gas dynamics, it 
is possible to justify it through the hydrodynamic limit (see a review on these questions 
e.g. in [23]). Here, the hydrodynamic limit can be developed solely in the special case 
where the interaction is local (as in section 13.51) and in the absence of any blind zone 
behind the subject. We will first propose two other closure methodologies which apply to 
general cases, but which cannot be justified by a hydrodynamic procedure. 

The first closure scheme, referred to as the 'monokinetic closure', is developed in sec- 
tion U]2j It is valid when there is rigorously no noise (i.e. no uncertainty in the pedestrian 
velocities). It postulates a monokinetic distribution function: in the neighborhood of a 
given location x at time t, all pedestrians having the same target direction a have the 
same velocity U(x,a,t). In other words, in this neighborhood, the statistics of possible 
velocities is given by a Dirac delta in the velocity variable, located at U(x,a,t). The 
resulting CM belongs to the class of second-order models of traffic: it involves two bal- 
ance equations for the mass and momentum densities respectively and bears analogies 
with pressureless gas dynamics models [13l H3]- These models have somehow unpleasant 
features, such as the possible formation of mass concentrations. 

For this reason, a second closure scheme, referred to as the 'VMF closure', is proposed 
in section H~3l The model supposes that some noise is involved in the pedestrian velocities; 
this is indeed more realistic than the zero-noise assumption of the previous closure. The 
distribution of velocities is supposed to be a von Mises-Fisher (VMF) distribution. The 
VMF distribution is a natural extension of the standard Gaussian distribution for random 
variables defined on the sphere [57]. Like in the monokinetic closure scheme, the resulting 
CM belongs to the class of second-order models of traffic but the form of the momentum 
equation has not been previously found anywhere else. 

Finally, in section |4.4[ we develop the hydrodynamic limit in the special case of a local 
interaction with no blind zone. The hydrodynamic limit supposes that, for a pedestrian, 
the process of turning towards the velocity which minimizes the potential $ is very short. 
Therefore, the velocity distribution can be approximated by an equilibrium which reflects 
the instantaneous equilibrium between the turning process and the noise. Such a distribu- 
tion, which will be our closure Ansatz, is called a 'Local Thermodynamical Equilibrium' 
(LTE), by analogy to the standard terminology of statistical mechanics. The LTE is very 
peaked around the velocity which minimizes the potential, with some spread due to the 
noise. An important point is that, while the LTE depends on the potential, the potential 
also depends on the LTE through the definition of the DTI. Therefore, the allowed DTI 
are determined by a fixed point equation. The resulting LTE can be interpreted as a Nash 
equilibrium of a game consisting for the pedestrians in finding the best compromise be- 
tween reaching their target and avoiding collisions with other pedestrians. The framework 
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for a game-theoretic interpretation of LTE can be found in [28]. The resulting model is a 
first-order model, in the traffic terminology sense. It consists of a conservation equation 
for the mass density, while the mass flux is determined functionally from the LTE, i.e. 
from the DTI that have been found by solving the consistency equation. 

4.2 Monokinetic closure 

4.2.1 Monokinetic closure: derivation 

In this section, in order to derive a macroscopic model, we assume a monokinetic distri- 
bution function. For the monokinetic assumption to be valid, we need to remove the noise 
term, and consider the following equation: 

d t f + cu ■ VJ + V u • (Ff) = 0, (4.3) 

coupled to (13.1 Op . The monokinetic closure consists of the Ansatz: 

f(x,u,a,t) = p(x,a,t)6u( x>ait) (u), (4.4) 

where Su(u) is the Dirac delta located at U (see a graphical representation at Fig. [8] (red 
arrow)). Note that, by definition, U(x,a,t) G S 1 i.e. is a vector of norm 1. This Ansatz 
means that there is only one definite velocity U(x, a, t) at any given point x, time t, for 
any preferred direction a. It is easily shown [23] that the distribution (j4.4p is an exact 
solution of (14 .3p provided that p and U satisfy the following set of macroscopic equations: 

dtp + V, • (cpU) = 0, (4.5) 

d t U + cU-V x U = F(x,a,t), (4.6) 

with F(x, a, t) = F(x, U(x, a, t), a, t) and F given by (13.101) . In other words, 

F(x,a,t) = -V w ® {Xta , t) (U(x,a,t)), with $( x , a ,t)(w) = ®(x,u(x,a,t),a,t){w), (4.7) 

and $ given by (13. lip . The potential $>( x ,a,t)(w) can be written: 

k 

®(x,a,t)(w) = -\D {x ^ t) {w)w-La\ 2 , (4.8) 

where D^ a ^{w) = D(x,U(x,a,t),w,t) is given by: 

__! ,, (SyeS(x,a,t)Sb&D~ 1 (v- x MyM-w)p(y,b,t)dvM i\, An . 

D, Aw) = max < — — T . T . — - — r— — , — ■ > (4.9) 

(x ' a ' t] 1 Iyes( x , a ,t)l^p(y^t)dydb 'L/ 1 



with 



S(x, a, t) = {y e M 2 | \y - x\ < 6(x, t), y- — ^- ■ U(x, a, t) > «}, (4.10) 

\y-x\ 



and the functions D(y — x,v — w) and 8(x, t) still given by (13. 14}) and ( 13. 17(1 . 

We note that, by definition (14.71) . F is orthogonal to U. Then, multiplying (14.6(1 
scalarly by U and using that the operator d t + cU ■ V x is a derivative, we get: 

(d t + cU-V x )(\U\ 2 ) = 0. 

Therefore, the constraint \U\ = 1 is satisfied at any time provided it is satisfied at time 
t = 0. 
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Direction of U 
Dirac delta distribution 6u{u) 



VMF distribution M v {u) 




Figure 8: The VMF distribution (in blue) and the Dirac delta distribution (in red) as 
functions of u in polar coordinates. The direction of the mean velocity U is given by 
the black semi-line. The width of the VMF distribution about its maximum u = U is a 
function of /3(\U\). In both cases, the mean velocity U is a function of (x, a, t) determined 
by the fluid model. We have \U\ = 1 in the Dirac delta distribution case and \U\ < 1 in 
the VMF case. 



4.2.2 Monokinetic closure: discussion 



The model expresses the conservation of mass (J4.5P and evolution of velocity (j4.6p . The 
mass conservation equation (14. 5 p takes the form of a classical continuity equation and 
expresses that the rate of change of p in any arbitrary small volume is solely due to the 
mass flow across the boundary of this volume. The velocity equation (14. 6p expresses that 
the rate of change of U along the flow lines (i.e. the left-hand side of ( 14. 6 p which takes 
the form of a material derivative of U) is proportional to the force F exerted on these 
particles. The target direction a does not explicitly appear in (I4.5| ) and ( I4.6P except 
through this force term which couples all target directions altogether. 

The force term (14. 7p describes how the bulk fluid velocity U changes in time: it 
tends to decrease the potential (14. 8p . whose minima express the best satisfaction of the 
target direction while avoiding collisions. The potential is computed as follows. By the 
monokinetic assumption ( I4.4p . all pedestrians within a given fluid element which have 
the same target direction a also have the same velocity U(x,a,t). Then, the elementary 
DTI of these pedestrians with pedestrians located at y and having target direction b is 
computed. Again, by the monokinetic assumption (14. 4p . these particles have velocities 
U(y,b,t). Therefore, this elementary DTI computed with the test velocity update w is 
given by D(y — x, U(y, b,t) — w), where D(y — x,v — w) is the elementary DTI (see 13. 14]) . 
Then, these elementary DTI are averaged over particles located in the cone of vision of 
x, defined at (j4.10p . This provides the averaged DTI D^ x>a> ^(w) with test velocity update 
w (formula 14. 9p . The average DTI is used to construct the potential ^( x , a ,t){w) which 
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expresses how far the pedestrian is from his target point when walking in the direction w. 
The expression of the force term is non-local in space: it involves the complex average 
(14. 9 p over a neighborhood of the point where the force evaluation is made. This non- 
locality expresses the ability of the pedestrian to anticipate the likelihood of a collision 
with the neighbors. However, a local version of this model can be designed, based on the 
local version of the kinetic model of section 13.51 This local version is simply obtained by 
replacing (14.91) by its local version issued from (13.241) . It leads to 

( f beSl A-i(U{x,a,t),U{x,b,t)-w)p(x,b,t)db 

D, Aw) = max < '^—^ ? ; — - — , 

(x ' M)l ' \ f b€S1 p(x,b,t)db 

with the function A Kt g(u, v — w) defined by (I3.23p . As seen in section |3T5| in the special 
case k = — 1 (i.e. there is no blind zone behind the subject), the function A Kt s(u,v — w) 
is replaced by A^jf — w\), which does not depend on u. In this case, the DTI does not 
depend on a and is given by the expression: 

n _! , , j/ks'^mG^m) -w\) P (x,b,t)dh n 

D « w = max 1 u**m«> • A • ( ] 

In all these cases, the evaluation of the force still requires an integration in the target 
direction variable. 

Apart from the complex expression of the force term, this model belongs to the class 
of pressureless gas dynamics models [131 El- Such models have some pathologies: they 
do not guarantee that the monokinetic closure assumption (see Fig. [H] (red arrow)) is 
preserved in time. Specifically, particle trajectories with same target direction a but 
initially located at different positions Xq and x' can meet at later times. This results 
in the appearance of a discontinuity of U (because the two meeting particle trajectories 
may have different velocities) and the blow-up of p. This classical phenomenon is similar 
to the appearance of caustics in geometrical optics. The non-local force term F(x, a, t) 
at the right-hand side of (14. 6 p is likely to be too weak to repel the trajectories at close 
encounters. To prevent this blow-up, it is necessary to introduce some kind of 'internal 
energy'. This is the motivation for the VMF closure below. 

4.3 VMF closure 

4.3.1 VMF closure: derivation 

In the previous section, it was not needed to take the moments of the noiseless kinetic 
equation (14.31) . since the monokinetic Ansatz (I4.4p provides an exact solution of it. Here, 
we will consider an Ansatz which is a priori not a solution of the kinetic equation ( 13. 91) 
(but which hopefully is close to one). But as a counterpart, we will be able to take into 
account the noise term. With this aim, we take the first two moments of ( 13.91) . 

First, integrating (13. 9p with respect to u leads to the mass conservation equation in 
the same form as previously: 

dtp + V x ■ (cpU) = 0. (4.13) 
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Indeed, the Fokker-Planck equation (I3.9P is of the form 

dtf + V, ■ (cuf) = V u ■ (A 1 + A 2 ), (4.14) 

where 

A 1 (u) = -Ff, A 2 (u) = dVJ\ (4.15) 

are tangent vector fields to S 1 . Therefore, thanks to Stokes's formula, 

/ V u • {A x + A 2 ) du = 0. 

Jugs 1 

Since u, x and t are independent variables, integration with respect to u commutes with 
derivation with respect to t and x and, with the definitions (14.11) . (14. 2ft . the integral of 
the left hand-side of (J4TT4"|) leads to the left-hand side of (fl~T3D . 

We now turn to the equation for pU . Multiplying (I4.14p by u and integrating with 
respect to u leads to 

dt(pU) + V, ■ (c£) = f V U -(A ± + A 2 ) u du, (4.16) 

with the 2x2 tensor S defined by 

£ = / fu <g) udu, 

Ju&S 1 

and m ® m is a matrix of components (u <8> u)ij = UiUj. Thanks to Stokes's formula, we 
have: 



V„- (A 1 + A 2 )udu = - / (Ax + A^du. (4.17) 

U& 1 Jues 1 

In particular, for A 2 , applying Stokes's formula once more, we get 

(V u • A 2 ) udu = —d / V u / du = — rf / f udu = —dpU. 



The other integrals, namely £ and that related to A\ cannot be expressed analytically 
from p and U. In particular, £ involves second order moments of / with respect to u. To 
proceed, we need a closure assumption. 

By analogy with gas dynamics, we assume that the velocity distribution is a von Mises- 
Fisher (VMF) distribution about the mean direction of U. The VMF distribution is what 
generalizes the Gaussian measure to the circle and more generally, to spheres [67]. In the 
present context, its expression is given by 

M C7 (n) = iexp{/3( M -fi)}, ^ = 7^7, (4-18) 

The quantity Z is a normalizing constant such that M\j is a probability density on S 1 : 

Mudu=l. (4.19) 

■.es 1 
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Then we have 

exp{/3 (u ■ Q)} du. 



1& 1 
Introducing the angle 9 = (Q,u), we can write: 

I-2TX 

Z=l exp{(3 cos9}d9 = 2irI (/3). (4.20) 

Jo 

We recall that Ik{x) denotes the modified Bessel function of the first kind: 

i r 

4(x) = - / exp{x cos9} cos(k9) d9, Vx e R, VA; G N. 
7T Jo 

The constant /3 plays the role of an inverse temperature: If (3 is large, Mjj(u) is extremely 
peaked in the direction u = Q while if /3 is small, My(u) is almost isotropic. It will be 
determined later on. A graphical representation of the VMF distribution My(u) as a 
function of u in polar coordinates is given at Fig. |8] (blue curve). The function Mu(u) is 
maximal at u — U. Its graphical representation is bell-shaped in a neighborhood of this 
maximum with a width roughly proportional to y^5. 

Now, we assume the VMF closure Ansatz, namely that / is proportional to My and 
is written: 

f(x, u, a, t) = p(x, a, t)Mu( Xt a t t){u), (4.21) 

where p(x, a, t) and U(x, a, t) are the moments (14. ip and (14.21) of /. That ( I4.2ip satisfies 
T} is obvious in view of the normalization condition (j4.19p . However, that it satisfies 
^2|) requires a relation between /3 and \U\ as we see now. Eq. (14. 2p requires that 

- exp{fi (u ■ £1)} u du = U. (4.22) 

We decompose u onto the direction spanned by Q and its orthogonal: 

u = (u ■ n)n + (u ■ n x )n x . 

Inserting this decomposition into (j4.22p and noting that the first term is an even function 
of 9 and the second one, an odd function of 9, we find that (j4.22p is equivalent to: 

/ — exp{0 (u • SI)} (u • £1) du = \U\. 

or equivalently, to: 

1 f 2n 

exp{/3 cos9} cos 9 d9 = \U\. 



\U\. (4.23) 



Z., 



This equation can be put in the form 

hiP) 



W) 
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The left-hand side is a monotonically increasing function of /3 G [0, oo) onto [0, 1) [25] . 
Therefore, as long as \U\ < 1, there exists a unique such that (14.231) holds. We will 
denote this root by (3 = j3(\U\). 

That \U\ should be less than 1 is a consequence of the kinetic model. Since the 
microscopic velocities u are normalized, the norm of their average over the probability 
M\j should be less than 1. The parameter \U\ is the classical order parameter of nematic 
materials (see applications to swarm dynamics in [33| |25| [65]). If \U\ is close to 0, then, 
(3 is close to and the distribution is almost isotropic. This indicates a disordered state, 
with microscopic velocities u pointing in all possible directions, so that their average 
almost cancels out. On the other hand, if \U\ is close to 1, then [3 is very large and the 
distribution is strongly peaked about the mean velocity Q. This is the case where almost 
all microscopic velocities are pointing in the direction of Q. 

Since Z depends on j3 through (I4.2UI) and (3 depends on \U\ through f)4.23p . we will 
now write [3 = /3(|C/|), Z = Z(\U\). The VMF distribution (14.211) is now written (omitting 
the dependences of p and U upon (x, a, t) for clarity): 

/(«) = pMu(u) = jjj^^P {^ (« " U)} ■ (4-24) 

Now, with (I4.24p . the tensor S can be computed. Introducing again the angle 9, we 
have: 

fu®udu = — I exp{/3 {u ■ il)} (u <g> u) du 
les 1 % Jugs 1 



p f n ( cos 2 9 cos 9 sin 9 \ rn „, ,„ 

Zj (.cos^sin^ sin 2 9 )^v{^os 9} d9 




l+cos 26* 







HP)- 



2 




q 1-COS26 ) ex P{/3 cos 9} d9 



In the third line, the off-diagonal terms have been canceled out because sin# is an odd 
function of 9. Therefore, introducing 

7„(CT~(l + m), 7i (|.7|) = ^(l -|gj), (4.25) 

with (3 = (3(\U\), we can write E as: 

£ = p(7||(|tf|) U ® U + 7±(l^l) ^ ± ®f/ ± )- (4.26) 

Since I2/I0 < 1> the matrix S is positive definite. The limit (3 — > is undefined since 
\U\ — >• and t^t has no definite limit. In the limit (3 — > 00, the matrix converges to 
pU (8> [/, and we recover the corresponding term of the monokinetic closure (second term 
at the left-hand side of (jOjl ). 
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The last term to be computed is the contribution of the force, i.e. the A\ term at 
the right-hand side of ( I4.17p . Using the first eq. ( I4.15P and (I3.10p . this contribution is 
written: 



Ai(x, u, a, t) du = p(x, a, t) F(x, a, t), 



u& 



F(x, a,t) = / F(x, u, a, t) M v ^ a>t ) (u) du, (4.27) 

with F(x,u,a,t) given by (j3.10p . The quantity F(x,a,t), which has a similar physical 
interpretation as in the monokinetic closure case (see (14. 7ft ). is the total force exerted 
on a ensemble of particles at position x and time t having same target direction a. The 
computation of F(x, u, a, t) follows that of section [3731 and is not repeated here. The only 
change brought by the VMF closure is in (13.131) . which can be written: 

n _! , v fIyes(x,u > t)Ibes^ E ' 1 (y- x ^ w ^ u (y^ b ^))p(y^ b ^) d y db il , Ano , 

D <XUt)( W ) = maX \ f f 1 I 4-\ -J JU ' T f ' ( 4 - 28 ) 

with 

E- 1 (y-x,w,U)= D~ 1 (y-x,v-w)M u (v)dv. (4.29) 

JveS 1 

The quantity E~ x {y — x, w, U) can be computed once for all. 

We now summarize the macroscopic model based on the VMF closure. It consists of 
the two equations for the mass and momentum: 

dtp + V. • (cpU) = 0. (4.30) 

d t ( P U) + V x - (cp(j\\U^U + 7±^®^)) =pF-dpU, (4.31) 

together with the expression ( 14. 27ft of the force F. We have omitted the dependences of the 
coefficients 7 upon \U\. This is a system for p(x, a, t) and U(x, a, t). The target direction a 
appears implicitly through the expression of the force F which couples all target directions 
together. The properties of this system, and particularly, its hyperbolicity, will be studied 
in future work. 

4.3.2 VMF closure: discussion 

The interpretation of the two equations ( 14.301) . (14.31)) is the same as for the monokinetic 
closure (see discussion at section 14.2. 2D . Eq. ( I4.30J) expresses the mass conservation, while 
eq. ( I4.3ip describes how the mean velocity evolves in time. Compared to the monokinetic 
closure, Eq. ( I4.3ip is more naturally written in terms of the momentum pU and presents 
three major differences. 

The first one is the expression of the transport term V x - [cp (7y U®U + 7j_ U l ®U ± )') ) 
which, compared to (14. 6p . involves two terms. The first one, proportional to U <g) U is 
similar to the term involved in (14. pp . but is multiplied by a coefficient 7||(|t^|) which is less 
than 1. The second term, proportional to the tensor U 1 - <S> U L , is unusual in fluid models. 
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It arises from the VMF closure, itself justified by the fact that microscopic velocities are 
constrained to be of norm 1. For such fluids, significant differences from classical fluid 
dynamics models have already been found [21)1 |2T?] . Similar unconventional models have 
been found even if the microscopic velocities are not constrained to be of norm 1 when 
the particles are subject to a self-propulsion force [6]. The impact of this new term on 
the hyperbolicity of the model will be studied in future work. 

The second difference is the presence of a velocity damping term (second term at the 
right-hand side of (I4.3ip ). This damping term is produced by the noise. At the kinetic 
level, the dynamics induced by the velocity diffusion operator is that of the heat equation 
on the circle S 1 . It makes the velocity distribution function more and more isotropic. 
Under this dynamics, the average velocity tends to zero which is what the damping term 
at the right-hand side of (j4.3ip expresses. 

The third difference is in the computation of the force term F. It is constructed by 
averaging elementary force terms over the probability distribution Mu(u) (see ( I4.27P ). 
The potential of each elementary force term depends on the averaged DTI of the particles 
belonging to the corresponding fluid element. This averaged DTI involves the average of 
the elementary DTI over the probability distribution Mu(u) again (formulas (I4.28P and 
( I4.29P ). The result of this procedure is a non-local expression of the force involving both 
the mass and velocity distributions p and U. 

We can derive a local version of the interaction force, using the local kinetic framework 
of section 13.51 In this framework, the averaged DTI is given by: 

°<-> W = maX J U/fr.M)* ' l) ' (4 ' 32) 

with 

£-J(w, w,U)= [ A-i(u, v - w) M v (v) dv. (4.33) 

The function £ K j can be computed numerically. 

In the special case k — — 1 (no restriction of the interaction region to a forward vision 
cone), the function A K( 5 becomes independent of u and the formulas simplify into 

n _i , x f/ftesi^it)^'^' 6 '*))^' 6 '*)* l\ f A o A \ 

D,,Aw) = max < — - r — - — — - , — > , (4.34) 



E^(w,U)= / AjWv-wfiMuWdv, (4.35) 



with 



and A$(\v — w\) given by (I3.23p . Again, the function 8j 1 {w 1 U) can be numerically 
computed once for all. We note that D7} t Jw) does not explicitly depend on u. It results 
that the potential does not depend explicitly on u either, and can be written 

k 
§{x,a,t)(w) = -\D {X)t) (w)w - La\ 2 . (4.36) 
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The elementary force can now be expressed by a gradient with respect to the u variable: 

F(x, u, a, t) = -V u $ (Sj0)t) («). (4.37) 

This was not possible before. Indeed, because of the explicit dependence of $ on u, we 
had to distinguish between the variables u and w in the expression of the force (I3.10p . 
In the present case, thanks to (I4.37p . eq. (I4.27P for the force F(x,a,t) can be simplified. 
Indeed, using Stokes formula in (I4.27P together with (I4.37p . we get, for each component 
k = 1, 2 and denoting by e k the k-th basis vector: 

F k (x,a,t) = -(/ V u ®{x,a,t){u)M U ( x>ait ){u)du)e k 
\Jues 1 J 

^u^(x,a,t)(u) M U{X:a>t )(u) [e k ■ U^)u^ du 

ues 1 
&(x,a,t)(u) V u • (Mufaafiiy) u^u 1 -) du, 



F(x,a,t) = P(\U{x,a,t)\) / ®( x ,a,t)(u) M u{x ^ t) (u) du 



which, after easy computations, leads to 

U(x, a,t) 
'ues 1 \U{x,a,t)\ 

U(x, a,t) 
'ues 1 \U(x,a,t)[ 

Inserting (14.361) into (14.381) leads to an alternate expression of the force (in the case 

K= -1). 



$( Xl a,t)(u) M u{x ^ t) (u) (l +P(\U(x, a,t)\) (u ■ i^'°'* , )) udu. (4.38) 



4.4 Hydrodynamic limit 

4.4.1 Hydrodynamic limit: derivation 

Here, we consider the situation of section 13.51 We assume that the various length scales 
associated to the interactions between pedestrians are very small. Additionally, we assume 
that k = — 1, i.e. the interaction region is the disk B$ and there is no blind zone behind 
the subjects (see (|3.25jl ). In this case (see bottom of section 1331) . the potential does not 
explicitly depend on u. In the present section, we take advantage of this simplification 
to perform the hydrodynamic limit of the mean-field model (13. 9p . The hydrodynamic 
limit consists in supposing that the changes in pedestrian velocities due either to the 
interaction force F or to the noise with diffusion constant d occur on very short time 
scales. To highlight this scaling assumption, we introduce a small parameter e< 1 and 
rescale the force and diffusion constants in such a way that 

F = -F, d = -d. (4.39) 

Dropping the 'hats' for simplicity and writing all unknowns with a superscript e, we 
can write the mean-field model (13.91) with local interaction force and no blind zone as 
follows: 

d t f + cu-V x f = -Q Dfe (r), (4.40) 

e 
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where the operator Qo f (f) (the so-called 'collision operator' of kinetic theory) describes 
the rate of change of the pedestrian velocities due to their interaction with the other 
pedestrians and the noise. We parametrize the collision operator by the DTI Df. Be- 
cause of the local interaction assumption, Qof{f) operates only on u and a, leaving (x, t) 
untouched. Therefore, we describe it as an operator acting on functions of f(u, a) only. 
For a given function u G S 1 — > D(u) G R+, the collision operator is written 

QdU) = -V„ • (F D f) + dAJ, (4.41) 

where the force Fd(u, a) is written in terms of the potential $d(w, a) by: 

k 
F D (u,a) = -V u $ D (u,a), $ D (u,a) = -\D{u)u - La\ 2 . (4.42) 

For a given function f(u, a), Df(u) is defined by 

, S {v , b) ^^(\v ~u\) f(v,b)dvdb 

D f \u) = max ^ , T } , 4.43 

S{v,b)eT*fv°Mdvdb 



5 f = c( [ f(v, b) dv db) 1 2 . (4.44) 



and As(\v — u\) is the known function given by (13.261) . Finally, Sf is given by 

-1/2 

f(v, b) dvdb] 

'(v,b)£T 2 

It is an immediate matter to check that this sequence of definitions is equivalent to the 
mean- field model ( 13. 9ft with local interaction force (13. lOf) . ( 13. lip and ( 13.24p . in the case 
K = — 1 (up to the change of scale ( 14.391) ). Here, we note that the first equation ( 14. 42ft 
is equivalent to ( I3.10p . because in the case k = —1, the potential & x ,a,t( w ) d° es no ^ 
explicitly depend on u. These two formulas are not equivalent if the potential (13. lip 
depends explicitly on u, which is the case when k > — 1. This is why this section is 
restricted to the case k — — 1. 

Now, by inspecting (I4.40H in the limit e — > 0, we deduce that f e — > f° where /° is the 
solution of 

QD /O (f°) = 0. (4.45) 

Any such solution is called a Local Thermodynamical Equilibrium (LTE). By the fact that 
the operator / — > Qo f (f) only operates on (u,a), this equation specifies the dependence 
of f° on (u,a), leaving the dependence on (x,t) undetermined at this level. 

In order to determine the LTE's, we first suppose that D: u G S 1 — > D(u) G IR+ is a 
given function. We note that for a fixed function D, the operator / — > Qn{f) is a linear 
Fokker-Planck operator. We introduce the function 

M D (u,a) = — — exp ( - D \> } ), (4.46) 

with Zo{a) the normalizing constant such that 

M D (u,a)du=l. (4.47) 



«e 



;i 
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The function m6§'4 Md(u, a) for a given a 6 S 1 is represented graphically on Fig. [9] 
(blue curve) as a function of u in polar coordinates. We notice that Mp is maximal where 
the potential $£> (represented by the black dashed curve) is minimal. Around its maxima, 
the graphical representation of Mp is bell-shaped. The corresponding widths are roughly 
proportional to the noise level yd. 



Direction of U , 



LTE distribution M D {u) 




Direction of the global 
minimum of & o 



Potential $>£> («) 



Figure 9: The LTE distribution u G S 1 — > Me>(u, a) for a given a G S 1 as a function of u in 
polar coordinates (blue curve). The distribution Mo is maximal where the potential $£> 
(black dashed curve) is minimal. The minima of $d and maxima of Mp are indicated by 
black semi-lines. The global maximum of Mp corresponds to the global minimum of $d- 
The distribution Mp has bell-like shapes around its maxima. Their width are roughly 
proportional to the noise level yd. The direction of the mean velocity U is given by the 
red semi-line. It is fully determined by Mo and therefore, by &d and is a function of 
(x,a,t) like $£>. We have \U\ < 1. 



Thanks to (14.42ft . we can write: 



QdU) = -dv u - {m d v u {^-)) 



(4.48) 



By applying Green's formula, we deduce that for any function f(u,a) with appropriate 
regularity, we have: 



QD(f)ji-duda = - 

M D 



M 



D 



(w,a)GT 2 



V " ( ^> 



duda. 



(4.49) 



'(u,a)£T 2 

Therefore, the equation 

QdU) = 0, (4-50) 

is equivalent to saying that there exists a function p: aeS 1 -) p(a) > such that 

f{u,a) = p{a)M D {u,a). (4.51) 
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Indeed, inserting (j4.5ip into (J4.48P clearly leads to ( I4.50p . Reciprocally, if ( I4.50P is satis- 
fied, then, applying (J4.49P leads to the fact that V U (//M D ) = 0, i.e. f/M D is a function 
of a only, which is exactly saying (14.5 ip . 

Therefore, an LTE is necessarily of the form (14.511) . However, there is a consistency 
condition to be satisfied. For / = pMo to be a solution of (14 .45 p . we need to ensure 
that D = Df. In other words, D is not arbitrary, but must be the DTI associated to /. 
Inserting (j4.5ip into (14. 43 p . this constraint is written: 

D («)=n«< j-^ , T >. (4.52) 




For any function p: a G S 1 — >• p(o) > we look for functions -D p : u G S 1 — > D p (u) > 
of this functional equation. The questions whether it admits a solution and how many 
such solutions exist is left to future work. Here, we assume that there exists at least one 
isolated branch of solutions D p . Therefore, the LTE are of the form pMp , with D p a 
solution of ( I4.52p . By restoring the dependence upon (x, t), we conclude that the solutions 
of (I4.45P are of the form 

/°(x, u, a, t) = P( x ,t)(a) M D (u, a), (4.53) 

where, for any (x,t), the function u — > D p (u) satisfies (14.521) with p(b) replaced by 
P(x,t)(b). By the normalization condition (I4.47|) . p(x,t){ a ) appears as the density of pedes- 
trians at point x and time t with target velocity a. It has the same meaning as p(x, t, a) 
in the monokinetic or VMF closures, but is written differently to highlight its dependence 
on the target velocity a. Indeed, the LTE at point (x, t) depends functionally on the 
function a — > p( Xjt )(a). By contrast, the LTE at (x,t) does not depend on the function 
(x,t) — > p( Xtt )(a)- This motivates this dissymetric treatment of the dependences of p on 
a on the one hand and (x, t) on the other hand. Now, we are looking for the equations 
specifying how p( x ,t)(o) varies with (x,t). 

To do so, we first notice that, because of the divergence form of (14 .48 p . we have 



/ Q Df (f)du = 0. 



(4.54) 



Therefore, integrating (J4.40P with respect to u and using (j4.54p leads to 

d t p £ + V, • (cp e U £ ) = 0, (4.55) 

with p e (x, a, t) the density and U e (x, a, t) the mean velocity of pedestrians at position x, 
time t and target direction a, given by: 

p e (x,a,t)= / / e (x,w,a, t) du, (p e U £ )(x,a,t) = / f e (x,u,a,t)udu. (4.56) 

Eq. (I4.55P is the continuity equation for pedestrians of target velocity a and is valid all 
the time (i.e. even when e is not small). Now, taking the limit e — > in (I4.56P and using 
the fact that f £ — ¥ f°, where /° satisfies (I4.53p . we get 

p £ (x,a,t) ->P( Xj t)(o), U £ (x,a,t) ^U P(xt) (a), (4.57) 
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with 

U p {a)= / M Dp (u,a)udu. (4.58) 

Of course, taking the limit e — ¥ in the continuity eq. (I4.55P leads to 

d t p(x,t)(a) + V x ■ (cp( x ,t)(a)U P{x t) (o)) = 0. (4.59) 

To summarize, the hydro dynamic model provides the spatio-temporal evolution of the 
density p(x,t)(o) °f pedestrians with target direction a. It consists of the single continuity 
eq. (14.590 . The different target directions are coupled together through the computa- 
tion of the average velocity U p (a) by means of ( 14.58)) . It requires the determination of 
the DTI at the same point. The DTI D p {u) is the solution of the functional equation 
( 14.52)) . parametrized by the function p(a). This computation is local in space-time but 
must be realized at any discretization point in space-time (x,t). Therefore, the practical 
determination of the velocity U p t) (a) may require high computational power. However, 
the local character of the problem is perfectly adapted to massively parallel or graphical 
computers. Note that this hydrodynamic model belongs to the class of first order models 
of traffic, since the velocity is fully determined by the knowledge of the density. 

4.4.2 Hydrodynamic limit: discussion 

The rationale of this model is best understood if a time discretization is performed. 
Suppose that the distribution p x (a) of pedestrians at position x and target direction a 
is known at time t n = nAt. We update this density at time t n+l using the following 
time-discrete version of the continuity eq. (14.590 by the Euler method: 

P n x +1 (a) = p n x (a) - At V, ■ (cp n x (a)U pnx (a)) = 0. (4.60) 

To use this scheme, it is necessary to compute the velocity U p ^{a). For this purpose, 
the DTI D p n{u) need to be computed by solving the functional equation (14.52)) . where 
p n is substituted for p. Once D p n is known, the LTE (14.530 can be computed and as a 
by-product, the mean velocity U p n (a) through ( 14.580 (see Fig. |9j the mean velocity is 
represented by the red semi-line). 

The heart of the model is the process of finding the velocities, given the density p n 
of pedestrians having prescribed target velocities. This process is decomposed as follows. 
First, knowing the density p n , the DTI in all directions are computed by solving the 
functional equation ( 14.52)) . We note that the DTI D p n(u) is independent of the target 
direction a and only describes the ability of a pedestrian to move in the direction u, given 
the density p n . The functional equation (14.520 describes how each pedestrian optimizes his 
actual velocity, i.e. minimizes the potential &d p , taking into account all other pedestrians 
around. This functional equation expresses a Nash equilibrium of a non-cooperative 
anonymous game with a continuum of players. Such games are characterized by an infinite 
number of players forming a continuum jl]. They are non-cooperative i.e. they exclude 
the possibility for the players to cooperate to improve their gain [52] • Finally, they are 
anonymous in the sense that two players with the same strategy cannot be distinguished 
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[63] . Recently, this category of games has been at the heart of the theory of 'Mean- 
Field Games' [17]. Here, the strategy variable of the players is the velocity u, while the 
target direction a is the players' type (see an introduction to game theory in [H]). The 
players' utility function is the opposite of the potential 4>d(u, a). The functional eq. (j4.52p 
expresses that each pedestrian separately cannot improve his utility function by choosing 
a different velocity u, which is the definition of a Nash equilibrium. This model is a 
particular example of the framework relating game theory and kinetic theory developed 
in [2S]- This viewpoint will be further expanded in future work. 

We stress the local character of the model: this optimum is realized locally, i.e. at 
any point x and at all times t. Once the equilibrium D p n has been found, the LTE and 
the mean velocity U p n follow directly. We note that the dependence of U p upon p is 
functional, i.e. the value U p (a) for a given target velocity a depends on p(b) for all values 
of b e S 1 . This can be understood easily. If there are more pedestrians heading towards 
a given direction, say bo, the DTI will be affected in all directions u and correlatively, the 
mean-value of the velocity of pedestrians heading towards direction a will be changed, 
even if a is very different from bo . 

5 Comparison with previous work 

In this section, we compare our results to the literature. We refer to [521 [53] for a discus- 
sion of the original discrete IBM. The time-continuous IBM (section I3.2J) obviously bears 
analogies with the social force model [351 [371 [38] • However, in our model, the elementary 
binary interactions are combined in a non-linear way (i.e. they are nonlinearly additive, 
see [H1IS3] for a discussion of this point). The velocity potential (I3.6P is reminiscent of the 
'steering potential' model of [H]. However, in |4"T] . the potentials of the various obstacles 
are added linearly by contrast to the present work, as already mentioned. Analogies also 
exist with the optimal control model of [10] . Indeed, in our work, the potential is similar 
to a cost function that the pedestrian dynamic tends to minimize. In [ID], three types 
of costs are considered: (i) the cost of drifting away from the planned trajectory, (ii) 
the cost of walking too close to other pedestrians and (iii) the cost of acceleration. In 
our constant velocity model, we have not considered any cost associated to accelerating 
(i.e. turning), which is probably incorrect. We also replace cost (ii) by a constraint (the 
distance traveled in direction w cannot exceed the DTI D^w)), and we minimize cost 
(i) subject to this constraint. The smoother expression of cost (ii) in [ID] allows for the 
inclusion of advanced features, such as the body compressibility. However, adding more 
features increases the number of parameters that need to be calibrated from the data. 
Our model has quite few parameters which need to be calibrated, which is an advantage 
in the context of scarce noisy data. 

We now turn to the mean- field kinetic model (section I3.3p . As already mentioned, 
kinetic models are scarce in the literature [DJ. Ref. [7J proposes a general kinetic framework 
for traffic and crowd dynamics but the specific features of pedestrian interactions are not 
detailed. The model presented in section 13.31 seems to be one of the very first crowd 
kinetic models based on a detailed analysis of pedestrian behavior. 

Most fluid models for crowds have been envisioned as extensions of road traffic mod- 
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els. Fluid models for traffic roughly fall in two categories: (i) first-order models which 
are composed of the continuity equation and an algebraic equation relating the flux to 
the density and (ii) second-order models, where the continuity equation is complemented 
with an evolution equation for the mean velocity. The prototype of first-order models is 
the Lighthill-Whitham- Richards (LWR) model [19] , while second-order models are repre- 
sented by the Payne- Whitham (PW) [58] and Aw-Rascle (AR) [5] models. Clearly, our 
first two fluid models, the monokinetic closure (section I4.2[) and the VMF closure (section 
I4.3p belong to the class of second-order models, while the third one, obtained through the 
hydrodynamic limit (section I4.4J) . is a first order model. 

We first discuss the hydrodynamic model of section 14.41 in reference to the literature 
on first-order models. The difficulty with extending the LWR model of traffic to crowds 
is the passage from one to two dimensions in the prescription of the flux. While the 
traffic flux is a scalar quantity, the pedestrian flux is a vector and a prescription is needed 
to decide of its direction. In this respect, there are three classes of first order crowd 
models: (i) those where the direction of the flux is fixed locally as a function of the 
density or its gradients jTOj [19], (ii) those where the direction is fixed by a non-local 
average [60] and (iii) those where the direction is given through the solution of a Partial 
Differential Equation (PDE), such as the Eikonal equation [121 SSI HI]- Our model bears 
the strongest analogies with the third type. Indeed, in the determination of the flow 
direction the Eikonal equation is replaced by the functional equation (14 .43 p . which offers 
a similar level of implicitness. Additionally, in [43] , the pedestrians minimize their travel 
times through an estimate which gives more weight to crowded areas. Therefore, the 
optimization principles underlying the dynamics of both our model and [13] are similar. 
However, in [13] , only the total density is taken into account in the travel time estimation, 
while our model also includes information about the velocities and target velocities. In 
particular, the following behavior (i.e. the fact that a pedestrian does not necessarily see 
a pedestrian moving in the same direction as himself as an obstacle) is likely to be better 
taken into account in our model. 

Second-order models for crowds are more scarce. Basically, referring to the classifi- 
cation in the paragraph above, only classes (i) (local prescription of the direction of the 
flux within the PW model [1] or AR model [8]) and (iii) (coupling of a PW model with 
an eikonal equation for the flux direction [16]) have been explored. The models of [36|[39] 
may be seen as belonging to class (i), although they involve a third equation (the energy 
balance equation). Our monokinetic and VMF closure models (sections 14.21 and I4.3P bear 
the strongest analogies with these last two references as they are obtained using similar 
methodologies (by closing a moment hierarchy from a kinetic equation). The least differ- 
ence arises in the case of the monokinetic closure, which is close to a zero-temperature 
fluid equation. The non-local expression of the force or its local approximation still make 
the specificity of our model. In [3^1 EH] , the effects of the interactions between the pedes- 
trians are mostly embedded in the energy balance equation. The VMF closure offers 
more differences with the fluid mechanical view of [36| [39] . Indeed, the latter postulate a 
Maxwellian (i.e. Gaussian) velocity distribution of pedestrians. This is obviously unlikely 
since pedestrians cannot reach arbitrarily large velocities. Our microscopic dynamics, 
which constrains the velocities to be of constant norm, more closely mimics the real be- 
havior of a pedestrian (at least at moderate densities where he is able to walk most of the 
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time at his free speed). As a consequence of this velocity constraint, the closure cannot be 
Maxwellian, but is a VMF distribution instead. We have seen in section 14.3.21 that very 
specific features emerge from this unusual 'hydrodynamics'. Obviously, these features are 
not taken into account in standard fluid models. 

Finally, the last point of our discussion is a comparison between the three models 
presented in this paper. The third one, which is related to game-theoretic concepts, 
captures nicely the mechanisms by which pedestrians achieve a consensus and maintain 
traffic efficiency even in very crowded environments. However, it is restricted to local 
interactions and uniform vision (i.e. no blind zone). Although actual pedestrians do 
have a blind zone, the approximation involved in the "no-blind zone" model may find its 
justification as follows: pedestrians walking behind the subject are unlikely to significantly 
modify the value of his DTI. Indeed the threat of a collision of the subject with these 
pedestrians is weak, since both are walking with the same speed. Therefore, the DTI 
is identical, weather its computation includes all pedestrians in the neighbourhood of 
the subject or excludes those who are inside the blind zone. Consequently, there is little 
approximation involved in replacing the actual vision cone by the entire neighbourhood of 
the subject. Considering interactions within their local approximation makes the model 
only suitable to the large-scales, where the interaction region of the subjects is small 
compared to the size of the scene. Another drawback of the model is the complexity 
of solving the fixed point equation for the DTI everywhere in space-time. Unless a fast 
algorithm is found and massively computers are used, this can lead to overwhelming 
computer costs. 

For this reason the second model, which relies on the VMF closure constitutes a 
good compromise between physical accuracy and computational efficiency. Its usability 
however is subjected to its well-posedness, i.e. to its hyperbolicity, a property which 
still remains to be investigated. Finally, the first model, which relies on the monokinetic 
closure, is the simplest one. At least, it guarantees local-in-time well-posedness, i.e. the 
model has a unique solution until two pedestrian trajectories issued from initially different 
points meet. At such a meeting point, a mass concentration occurs, which is obviously 
unphysical. A way out this unpleasant feature would be to restore non constant speeds 
and to allow pedestrians to slow down in case of a close encounter, which they obviously 
do in practice. Consideration of non-constant speeds is the subject of future work. 

6 Conclusion 

In this article, we propose a hierarchy of macroscopic models derived from the heuris- 
tic behavioral Individual-Based Model of [53] and discuss them in view of the available 
literature. This IBM supposes that pedestrians optimize their trajectory towards their 
goal subject to the constraint of avoiding collisions with neighboring pedestrians. We 
first propose a novel kinetic model. In a second step, we derive three different fluid mod- 
els. The first two ones consist of balance equations for the density and mean velocity of 
the pedestrians. They differ by the proposed closure relations, based on a monokinetic 
Ansatz for the first one and on the VMF distribution for the second one. The third model, 
which has more restrictive assumptions, consists of a single mass conservation equation, 
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with mass flux functionally related to the density distribution. The functional relation 
expresses the realization of a Nash equilibrium where each pedestrian finds his optimal 
direction of motion towards his target in the midst of the other pedestrians and in the 
presence of noise. These models are the first available kinetic and fluid models derived 
from the heuristic behavioral Individual-Based Model of [53]. Future work will be devoted 
to the study of these various models, both from the theoretical and numerical viewpoints, 
and to their confrontation with the experimental data. 
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